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Abstract 

Properties of hybrid stars with a mixed phase composed of asymmetric nuclear matter and strange 
quark matter are studied. The quark phase is investigated by the quark quasiparticle model with a 
self-consistent thermodynamic and statistical treatment. We present the stability windows of the 
strange quark matter with respect to the interaction coupling constant versus the bag constant. 
We find that the appearance of the quark-hadron mixed phases is associated with the meta-stable 
or unstable regions of the pure quark matter parameters. The mass-radius relation of the hybrid 
star is dominated by the equation of state of quark matter rather than nuclear matter. Due to the 
appearance of mixed phase, the mass of hybrid star is reduced to 1.64 M with radius 10.6 km 
by comparison with neutron star. 

Keywords: Thermodynamic and statistical physics of quasiparticle model; Equation of state of 
strange quark matter; Phase diagram of hybrid stars 



1. INTRODUCTION 

The appearance of quark matter or hadron-quark mixed phase in the massive neutron stars is 
a hot topic in the study of compact objects. The baryon densities of the stars cover a larger range 
from very low densities in the outer part to the order of about ten times the saturation density in 
the inner core. To study the structure of compact stars, the key point is to find a reliable form of 
the equation of state (EOS) that determines the characteristic of the constituent matter [1]. Unfor- 
tunately, however, there is no single theory to cover the large density range with respect to quark 
degrees of freedom. In literature, the low-density phase can be described by quantum hadro- 
dynamics (QHD). At high densities, a new form of matter, called strange quark matter (SQM), 
might exist and be more stable than nuclear matter ( 56 Fe) yl |j, |4|, |5l |6J] . For more theoretical 
and experimental results, see 10 S Isl • A family of compact stars consistin g co mpletely of 
the deconfined mixture of u-, d-, s- quarks has been called "strange stars" 111 ll \l2i 1 1 3fl . If the 
hypothesis of stable strange quark matter is correct [3i), the possibility of phase transition exists 
in principle [ 14]. The compact star with a "quark matter core", either as a hadron-quark mixed 
phase or as a pure quark phase, are called "hybrid stars" (HyS) lfl5l[l6ll . Recently, it was argued 
that the interior core of a low-mass compact star could be dominated by the Color-Flavor-Locked 
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(CFL) quark matter II17L 1 1 811 - In contrast, it was remarked that the CFL phase would make the 



HyS unstable ifisill and the possible configuration of compact stars, such as the strange hadrons, 
hyperonic matter 1I20I1 and quark matter core, can soften the equation of states of neutron stars 
I 20l 21, 22]. However, Alford et al pointed out the "masquerade effect" |23] that the hybrid star 
has a mass-radius relation similar to that of the pure neutron star. Also, the isolated neutron star 
RX J 1856-375 with a large radius and/or mass [24] is a possible candidate of HyS, which implies 
a constraint on testing a rather stiff equation of state at high density rf25ll . Regardless of whether 
the quark matter is ruled out, it seemed that the soft equation of state were ruled out in the center 
of compact stars ll26Tl . 

To get a reliable equation of state in the microscopic calculation of interacting dense hadronic 
matter, Wiringa et.al added the three-body potential to the nucleon Hamiltonian and gave the light 
nuclei binding energy and nuclear matter saturation properties B27I1 . The Brueckner theory with 
three-body forces has been used recently to study the mixed phase of hadrons and quarks in 
compact stars [28]. In literature, there is another successful method called relativistic mean-field 
(RMF) theory. It is a powerful tool in describing various aspects of nuclear physics, such as the 
properties of nuclear matter, finite nuclei, and neutron stars, as well as the dynamics of heavy- 
ion collisions [29, 30]. Recently, the model has been extended to include the density-dependent 
meson-nucleon coupling constant in finite nuclei [31] and nuclear matter 125, 321. 

In the theoretical description of the deconfined quark matter, we can resort to the phenomeno- 
logical models constrained by experimental information. There are many successful works con- 
sidering the medium-effect of quark masses. One of them obtains the confinement by the density 
dependence of quark masses 113311 . Another one is the quark quasiparticle model, where the vac- 
uum energy density is not constant but density-dependent. For the medium dependence of the 
quasiparticle masses, it was derived at the zero-momentum limit of the dispersion relations from 
an effective quark propagator by resuming one-loop self-energy diagrams in the hard-dense-loop 
approximation (HDL) [34]. We have recently extended the model to include the important finite- 
size effect [35] and magnetized strange quark matter (MSQM) ll36ll . 

In this paper, we study the effect of the interaction coupling constants on the equation of 
mixed phase of nuclear and quark matter and the properties of hybrid stars. For the nuclear 
EOS, we adopt the relativistic nuclear field theory solved at the mean-field level, and especially 
the Baym-Pethick-Sutherland (BPS) model for densities below 5 x 10 14 g cm" 3 |37l S]. In 
describing quark matter, we employ the quark quasiparticle model instead of the conventional 
bag model [22]. The density-dependent bag function is obtained self-consistently rather than 
artificially. 

This paper is organized as follow. In Section [2] we briefly give a short introduction to the 
relativistic nonlinear mean-field model describing the nuclear matter. In Section [3] we introduce 
the treatment of SQM in the framework of the quark quasiparticle model and present the stability 
window of quark matter. In Section |4] we display the phase diagram of the mixed phase and 
discuss the chemical potential behavior. With the equations of state, we investigate the influence 
of coupling constant and bag constant on the mass-radius relation. The last section is a short 
summary. 



2. The nuclear EOS in the relativistic mean-field model 

The relativistic nuclear field theory is solved at the mean-field level. The in-medium inter- 
action of nucleons can be realized through the exchanges of the <x, u, and p mesons. The bulk 
matter is assumed to be electrically neutralized and in the lowest energy state, i.e. in general 
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/^-equilibrium. The influence of the temperature can be neglected 113911 . The Lagrange density for 
this model is written as ll40l 41 ] 

-Crmf = Ud^cr + mlp,ff)-U(cr) + V(c,)-j[n, v Q.^+R; v R a ^ + F, y F^] 



+tj/ \iyd M - m N + g aN cr - g^y^ - -gpNy 11 ^ ■ r a - eyfQeAft ) iff 
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where the nucleon field iff has a form of column vector for protons and neutrons, and the field 
tensors are given by 
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The non-linear potential functions are contained for the meson <x and to as 



1 2 2 ^ 3 1 4 

-m^cr + -g 3 a + -g 4 cr , 

1 2 1 2 

V{co) = -m^af + -c^aj^af) . 



(2) 
(3) 
(4) 



(5) 
(6) 



The parameters mj</, nig-, m p , m w are the masses of nucleon, scalar meson cr, isovector-vector 
meson p, and isoscalar- vector meson a> respectively. In principle, they can be fixed algebraically 
by the properties of bulk nuclear matter, such as, the binding per nucleon, the saturation baryon 
density, the effective mass of the nucleon at saturation, and the compression modulus. The 
isovector p field vanishes for symmetry nuclear matter. 

The scalar density and conserved baryon number density are expressed with the use of Fermi 
integrals 
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where the effective nucleon mass is defined by M* = M - ga-NO~u, with v p and v n being the Fermi 
momenta of protons and neutrons respectively. 

The energy density and pressure of nuclear matter can be obtained from the energy-momentum 
tensor. Including the contribution of nucleons and mesons, the total energy density e HP and the 
pressure are ifBiH 
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If the electron contribution is included, P e = /4/(127r 2 ). For symmetric nuclear matter, the 
number density of protons and neutrons are equal, and correspondingly the Fermi momenta are 
also equal, i.e., v p - v„. In this case, the isovector meson has no contribution. 

From the equation of motion for nucleons, the Fermi energy of nucleons, or equivalently, the 
chemical potential ^ can be expressed as 

UN = guN^O + ~gpNP03T3N + Vv 2 + M* 2 . (11) 

In this paper, we choose four typical sets of parameters (TM1, NL3, BKA20, and TW-99) 
ll^PPfll . They stand for different stiffness equations of state. The messes for a, a>, and p mesons 
are m ff = 509MeV, m w = 782MeV, and m p = 770MeV with nucleon mass M = 939MeV. For 
getting the EOS in the low-density, i.e., the crust equation of state of compact stars, we can resort 
to the BPS model B37I1 in calculations. 



3. Quark quasiparticle model for quark matter EOS 



At high density, nuclear matter is expected to undergo a phase transition to a deconfined 
phase. If the quark chemical potential exceeds the strange quark mass, the system can lower its 
Fermi energy by converting the down quark into strange quarks. Recently, we have developed 
the quark quasiparticle model in studying strangelets 113511 and MSQM |36]. In the quasiparticle 
model, the effective quark mass following from the hard-dense-loop (HDL) approximation of 
quark self-energy at zero temperature is expressed as [34] 451. 



m : — h 

2 



2 7 2 

4 + 6/r 2 



(12) 



where m, is the corresponding current mass of quarks, g = yjAna s denotes the strong interaction 
coupling constant. The effective quark mass m* increases with g, m, and the quark chemical 
potential ju,. In this paper, we treat g as a free parameter in the range of (0, 5). For light quarks 
(m and d quarks), we take the current mass to be zero, and Eq. (fTZt is reduced to the simple form 
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The quasiparticle contribution to the thermodynamic potential density is given as ll35ll . 

q. = -g j^{ln[l+ e -^-^ (14) 



where e, iP = ^jp 1 + m* 2 and T is the temperature. At zero temperature, the integration can be 
calculated out to give 
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where m, and /z; are, respectively, the particle masses and chemical potentials, d, is the degener- 
acy factor with di = 2(spin) x 3(color) = 6 for quarks and di = 2 for electrons. When the variable 
set (T, V, {yu,}) is chosen as the independent state variables, the thermodynamic potential is the 
characteristic function. In the quasiparticle model, the total thermodynamic potential density can 
be written as 

£2 = Yj ["ite. m *) + B M + fi o, (16) 

i 

Where BjQii) is the medium-dependent quantity determined by thermodynamic consistency. A 
little later, we will see that B,(/i,) is given by an indefinite integration. When it is expressed by a 
definite integration, an integration constant is needed. Therefore, Bo is from the sum of relevant 
integration constants, and we treat it as a free input parameter. In the quasiparticle model, the 
particle number density should be of the same form as that of a Fermi gas with the normal particle 
mass replaced by the effective quasiparticle mass, i.e., 



On the other hand, we know from the fundamental thermodynamics that 
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Equating the last equality in Eq. ( TT8l with the first equality in Eq. (1171 1. we immediately have 
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where the derivative of the thermodynamic potential density with respect to the quark effective 
mass can be analytically expressed by 
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Accordingly, the pressure pQ p , energy density e^ p , and baryon density p QP for SQM at zero 
temperature are written as, 
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Because the current mass of light quarks is nearly zero, one can have an analytical expression 
for light quarks by combining Eqs. (fT3l and (fT9l Il35[l34fl . giving 
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where a — g/( V67r) and B = y/ 1 - a(g) 2 are ^-dependent functions. 
For massive quarks, B*(p) has been observed as the result i34ll . 
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where the coefficient a„ can be related to the coupling constant through four polynomial function 



as in 113411 . which is very different from the gaussian parametrization of the density-dependent 
bag constant ll46ll . 

With the above quark mass formulas and thermodynamic treatment, one can get the properties 
of bulk quark matter. According to the Witten-Bodmer hypothesis, it is required that the energy 
per baryon of two-flavor quark matter is bigger than 930 MeV in order not to contradict with 
standard nuclear physics, but that of symmetric three-flavor quark matter is less than 930 MeV. 
For the co-existence of nuclear matter and quark matter in the process of phase transition, the 
allowed values of model parameters should be chosen in the metastable or unstable regions. In 
Fig. Q] we present the stability windows of SQM in the B^ 4 -g plane. The area below the dotted 
line is forbidden because the energy per baryon e QP /p QP of 2-flavor quark matter is less than 930 
MeV. The e QP lp QP of 3-flavor quark matter is smaller than 930 MeV below the solid line and is 
smaller than 939 MeV in the narrow area below the dashed line. In the top area the strange quark 
matter is unstable. 




Figure 1: Stability windows in the parameter space (B„ , g). Below the dotted line is the forbidden area where two-flavor 
quark matter is absolutely stable. Above the solid or dashed line, SQM is meta-stable or unstable, and SQM is absolutely 
stable between the the dotted and solid curves. 



6 



4. Hadron-quark mixed phase and hybrid stars 



4.1. Transition from nuclear matter to quark matter 

We suppose that the compact star composed of quark matter core and nuclear matter surface. 
The nuclear matter includes protons, neutrons and electrons. The quark matter consists of a 
mixture of quarks (u, d, and s) and electrons. In the interface of the two phases, there may be a 
mixed phase of nuclear and quark matter, which is important in understanding the properties of 
hybrid stars. To describe the structure of the mixed phase, we define the fraction^ occupied by 
quark matter by x = V® p / (V® p + V HP ). Then we have the total baryon number density p, energy 
density e, and electrical charge Q as, 

Pb = (1-X)P HP +XP QP , (26) 
e = (27) 

Q = d-x)Q HP +xQ QP , (28) 

where Q = according to the bulk charge neutrality requirement. The critical density p^ p for 

pure nuclear matter is determined by x - 0- With increasing density, the deconfinement phase 

transition takes place. Consequently, quark matter appears in the mixed phase. When x increases 

op 

up to 1, the quark phase dominates the system completely at the critical density ■ 

From the quark constituent of nucleons, we have the chemical potential relations through the 
linear combinations, 

fip = 2yu„ + jj d , ii„ = ii u + 2n d . (29) 

Incorporated with the /^-equilibrium condition [ij = fi s = p. u + jj. e , there are only two indepen- 
dent chemical potentials, e.g. p. u and fi e . The Gibbs condition for the phase equilibrium between 
nuclear matter and quark matter is that the system should be in thermal, chemical, and mechan- 
ical equilibrium. Therefore, in addition to the common zero temperature, we also have the two 
conditions 

^ = V QP , (30) 

P HP (^,yU„,yU e ) = P QP Qlu,Md,Hs,f^e), (31) 

where P HP and P® p denotes the pressure in the nuclear phase and in the quark phase, respectively. 
They are given by Eqs. (TTOb and (OTT i separately. The electron is uniformly distributed in the 
system and the pressure is common for both phases. Eq. (l30t means that the baryon chemical 
potential in both phases are the same. It can be shown that Eq. ( T30l is equivalnet to Eq. J29l lf28ll . 
For a given total baryon number density, one can obtain the two independent chemical potentials 
and^ by solving the equations d26l >, d28l ) and OTt . 

In Fig. 12 the system energy per baryon is displayed as a function of the density for g = 3, 
B l J 4 = 150 MeV, and TM1 parameter set. The mixed phase (shaded area) exists in the range of 
about 1 ~ 3 times the nuclear saturation density. In this density range, the energy per baryon 
of mixed phase (solid line) is lower than that of both the pure nuclear (dashed line) and the 
pure quark phase (dotted line). So the appearance of the mixed phase in neutron star matter is 
energetically favored for a proper parameter in meta-stable or unstable regions. This observation 
is consistent with that in Ref. 12811 . If one chooses the absolutely stable parameter-sets (g, B l J A ), 
no mixed phase will exist and the hybrid star is unstable and will collapse to a strange star. In the 
middle density range, the mixed phase starts at the nuclear critical density where x — and 
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Figure 2: The energy per baryon as a function of density. 
In the density range of shaded area, the energy per baryon 
(solid line) is lower than that of either quark matter (dotted 
line) or nuclear matter (dashed line), g = 3, Bq' 4 = 150 
MeV, and TM 1 set are adopted. 



Figure 3: The electron chemical potential fi? (left axis) 
and the fraction of SQM^ (right axis). Parameters are the 
same as for Fig. [2] 



ends at the quark critical density pf where x — 1 ■ The critical points are marked with full dots 
in Fig. |2 By comparing the different equation of state for nuclear matter, we find that the solid 
lines between the p^ p and pj? p is shorter for harder EOS of nuclear matter. It is also observed 
that the quark critical density pj? p varies only slightly with the RMF parameters. 

To provide a better understanding of the Glendenning hypothesis of global charge neutrality, 
we give the electron chemical potential and the quark fraction versus the total density in Fig. [3] In 
the shaded area of mixed phase, the increasing contribution of quark phase is a complementarity 
to the decrease of electrons. The electron chemical potential fi e (dashed line) decreases rapidly 
from the maximum value 85 MeV to 20 MeV. And finally it becomes very small in the pure 
quark phase where the quark fraction x is unity, where the three-flavor quark matter occupies the 
whole space. 

4.2. The structure of hybrid stars 

With the above equation of state, we now study the structure of compact stars. As usually 
done, we assume that the hybrid star is a spherically symmetric distribution of mass in hydrostatic 
equilibrium. The equilibrium configurations are obtained by solving the Tolman-Oppenheimer- 
Volkoff (TOV) equation for the pressure P(r), the energy density e(r) and the enclosed mass 
m(r): 

dP(f) = Gm(r)e(r) [1 + P(r)/e(r)][l + 4nr 3 P(r)/m(r)] 
dr r 2 1 - 2Gm(r)/r 
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where G — 6.707 X 10 45 MeV 2 is the gravitational constant, r is the distance from the center of 
the star. The subsidiary condition is 

dm(r) 2 



dr 



= Anr l e{r). (33) 



Giving the stellar radius R, which is defined by zero pressure at the stellar surface, the gravita- 
tional mass is given by 



Jo 



M(R) = 4n I r 2 e(r)dr. (34) 
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Figure 4: The mass-radius relation of neutron stars. The 
maximum mass is marked with an asterisk on each curve. 



Figure 5: The mass-radius relation of hybrid stars with 
different parameter sets of the nuclear matter. The max- 
imum mass, marked with a full dot on each cure, do not 
change significantly with the nuclear parameters. 



Before investigating the properties of the hybrid star, let's give the mass-radius relations for 
compact stars composed of pure nuclear matter. The TOV equation is solved with the nuclear 
matter EOS introduced in Sec. [2] We take the four parameter sets in the calculation: TM1, NL3, 
BKA20, TW-99. The maximum masses of neutron stars are located in 2 ~ 3 times the solar 
mass M with the corresponding radii in the range (12.3 ~ 14.1 km) in Fig. |4] For harder EOS 
of the parameter set TW-99, the maximum star mass approaches the value 2.86 M , which is 
bigger than the mass 2.06 M Q with the parameter set TM1. Different with the self -bound stars, 
the equation of state of nuclear matter at lower density is calculated in the BPS model. 

To investigate the influence of model parameter on the mass-radius relation of hybrid stars, 
we should consider different EOS of the nuclear matter and strange quark matter. Firstly, the 
parameter set (g = 3, B l J 4 = 150 MeV) is employed on the quark side, while different parameter 
sets on the nuclear side are adopted in the calculation. The corresponding mass-radius relations 
are plotted in Fig. where the maximum star mass are marked with full dots. Even though 
the nuclear parameter is adopted for different sets, the maximum mass is very close to the same 
value 1.64 M . Doing calculations with other choices of the parameter Bo and g, we find the 



result that the maximum mass does not be significantly influenced by the parameter uncertainty 
of the nuclear equation of state. 

Now we fixed the TM1 parameter set for the nuclear EOS while changing the values of g 
and Z?o on the quark side. In this case the maximum mass of hybrid stars is calculated in Fig. [6] 
The distinction between Fig. [5] and [6] can be understood from the comparison of EOS. In the 
upper panel of Fig. [7] we can know that the narrow range on the energy density axis is affected 
by nuclear parameter sets, while a larger range is greatly dominated by quark matter parameter 
sets. For a bigger value of the maximum of hybrid stars, the corresponding density range of the 
mixed phase is wider. In the parameter range, the maximum mass is in the range (1.5 ~ 2.O4)M 
and radius in the range (9.64 ~ 12.4) km, which comprises the pulsar PSR J1614-2230 with 
1.97 ± 0.04M o B47I1 . For a completely relation between the maximum mass of hybrid stars and 
the parameter set of SQM, we show the contour plots of the maximum mass in the panel with 
i?y 4 and g on the horizontal and vertical axis in Fig. [8] According to the stability window shown 
in Fig [1] the maximum mass of hybrid stars has been marked on each line in the allowed area. 




Figure 6: The mass-radius relation of hybrid stars with 
various parametrization of quark matter. The maximum 
masses are marked with full dots. 



Figure 7: The pressure versus the energy density for the 
equations of state of hybrid stars. Figure (a) and (b) are, 
respectively, for Fig.|5]and|6] 



5. Summary 

We have studied the hybrid stars with mixed phase of nuclear matter and SQM. In the outer 
layer, the nuclear matter is described by the relativistic mean-field theory and crust EOS by 
the BPS model. In the inner core, the quark matter is investigated by the quasiparticle model. 
According to the Witten-Bodmer hypothesis, we present the stability window of quark matter 
in the coupling constant g versus bag constant B l J 4 panel. The mixed phase exists in the range 
1 ~ 3 times the nuclear saturation density, which is energetically favorable by comparison with 
pure nuclear matter or SQM. We point that the mixed phase of hybrid stars is located in the 
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Figure 8: The contour plots of the maximum mass of hybrid stars in the panel (B^' 4 , g) . Based on the stability window 
in Fig.[T] the maximums of hybrid stars are shown on each line in the allowed area. 



meta-stable or unstable regions of pure quark matter. Otherwise the hybrid star will collapse to a 
strange star and no mixed phase could exist. 

We show that the maximum mass and radius of hybrid stars are controlled by the equation 
of state of quark matter |48] rather than that of the nuclear model. Or more exactly in the 
present framework, they are mainly controlled by the coupling constant g and bag constant Bo 
in the quasiparticle model. For a smaller coupling constant, the maximum mass is 1.54 M 
with radius about 9.64 km. The larger coupling constant will broaden the range of the mixed 
phase and increase the maximum mass of a hybrid star up to 2.0 M . So the quark matter with 
strong coupling can meet the constraint at high density in compact stars l25ll and corroborate the 
"masquerade effect" in the previous studies B23I1 . For even high density, the possible CFL matter, 
especially the magnetized CFL matter 1I49I1 . should be considered, which will be investigated in 
our future work. 
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